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Abstract—Gas condensate or natural gasoline is a valuable 
liquid hydrocarbon mixture which could be utilized as fuel or 
converted to different petroleum products like gasoline, light oil, 
jet fuel etc. It should be processed so that its water, salt, and 
acid contents meets the required standards for storage and 
transportation. Its vapor pressure, determined as Reid Vapor 
Pressure (RVP), must lie within a range that light components 
don’t evolve as a separate gas phase in transport pipelines or 
storage tanks. The optimum value is usually 10 psia in summer 
and 12 psia in winter. In the present study, we investigated 
condensate stabilization using two methods of multistage flash 
vaporization and distillation (fractionation) within two case 
studies. The results show that depending on the properties of 
feed to be stabilized, one or both of the two methods could be 
practical and/or economic. 


Index Terms—Gas condensate, HYSYS, ASPEN PLUS, 
optimization, simulation, stabilization. 


I. INTRODUCTION 

The process of increasing the quantity of intermediate (C; 
to Cs) and heavy (Cs.) components in gas condensate is 
called condensate stabilization. The main purpose of this 
process is to reduce the vapor pressure of the condensate 
liquids to prevent the production of vapor phase upon 
flashing the liquid to atmospheric storage tanks. On the other 
hand, the scope of this process is to separate very light 
hydrocarbon gases, methane and ethane in particular, from 
the heavier hydrocarbon components (C3,). Stabilized liquid 
has a vapor pressure specification, as the product will be 
injected into a pipeline or transport pressure vessel, which 
has definite pressure limitations [1]. It is usually 10 psia in 
summer and 12 psia in winter case. It should be noted that 
acidic contents (i.e., HS, CO», mercaptanes etc), glycol, free 
water, and possible salts must be also removed to have a gas 
condensate suitable for end-use applications. There are two 
methods for stabilizing gas condensates: flash vaporization 
and fractionation (distillation). Fractionation is more efficient 
and provides better control over vapor pressure for tanker 
shipment [2]. Non-refluxed distillation column is simpler but 
less efficient than refluxed one. Since it doesn’t require any 
external cooling source, it is particularly applicable to remote 
locations [2]. Benoy and Kale [3] investigated gas condensate 
stabilization through three stage flashing, non-refluxed, and 


Manuscript received October 10, 2012; revised November 24, 2012. 

The authors are with the Petroleum and Petrochemical Engineering 
School, Hakim Sabzevari University, Iran (Tel.: +98-571-4003513, fax: 
+98-571-4003648, e-mail: n.moghadam@hsu.ac.ir). 


DOI: 10.7763/IJCEA.2012.V3.244 


461 


refluxed distillation with and without hydrocarbon recovery 
from the flashed gas by compression. They found that heat 
duty in multistage flashing is much higher than distillation 
while liquid recovery was less. Nevertheless, liquid recovery 
was independent of the technique applied when intermediate 
hydrocarbons were recovered from the flashed gas. 

Herein, we undertook two case studies in which two-stage 
flash vaporization and distillation were investigated. The 
main objective was to adjust the RVP of product at 10 psia 
(good for summer) while keeping the duties of heat 
exchangers as low as possible and obtaining as much as 
stabilized condensate as possible. The process simulations 
were performed by means of ASPEN HYSYS and ASPEN 
PLUS software. 


II. PROCESS DESCRIPTION 

NGL-3100 gas plant is under construction near the oil 
production unit of Cheshmeh-khosh located at the southwest 
of Iran for optimal recovery of C2, product from the 
associated gases of North Dezful oil fields. It’s planned to 
produce 46,000 bbl.d' Cz,, 1180 bbl.d' Stabilized Sour 
Condensate, 6.12 MSCMD Sweet Gas, and 0.62 MSCMD 
Sour Gas. 

Shahid Hashemi-nejad (Khangiran) Gas Refinery is 
located at the northeast of Iran. The feed to the refinery 
comes from sweet gas reservoirs of Khangiran and sour gas 
reservoirs of Mozdoran. The treatment of 41.5 MSCMD of 
sour gas, dehydration of 7.4 MSCMD of sweet gas, 
distillation of 360000 L.d”' of gas condensates, and recovery 
of 2400 ton.d“sulfur are the most important operations 
undergone in this plant. The feed to the stabilization units of 
Khangiran and NGL 3100 plants are mixtures of filter 
separators, slug catcher, and regeneration gas separators of 
which the main specifications are listed in tables I and I. 
Peng Robinson equation of state was applied for 
thermodynamic calculations of both gas and liquid phases. 
Since there was no data available on the boiling point curves 
(including ASTM and TBP distillation curves), the 
assessment of C19, and C,, oil cuts specifications was carried 
out by using their respective specific gravities (Sp.Grs.) and 
molecular weights (MWs) [4]. 

The final temperature and pressure of the stabilized 
condensate are governed by the storage conditions which are 
35.3°C (95.5°F) and 2 bar in NGL3100 and 37.8°C (100°F) 
and 0.96 bar in Khangiran case. 
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TABLE I: FEED PROPERTIES OF KHANGIRAN GAS PROCESSING PLANT 


Comp. 10 No CoO, HS COS C, 
“ae 0 0.02 747 11.20 0 11.3 
Comp. G G CG, mC, kc; Cs 
— 0.97 0.59 0.91 127.0 «as 
Comp. Co Cc Cs Co Cio+ 

= 1.83 10.41 6.90 6.90 37.86 

a i 
74 739 293.6 T1904 1557 0715179 


TABLE II: FEED PROPERTIES OF NGL 3100 GAS PROCESSING PLANT 


STABILIZATION UNIT 
Comp. H,0 N> CO, HS Cos Cc, 
(mol%) 0.00 0.05 2.07 5.31 0 9.07 
Comp. C, C; 1-C4 n-C,4 i-Cs n-Cs; 
(mol%) 13.76 20.5 5.67 15.84 6.49 6.23 
Comp. Cos CS, a 
(mol%) 14.92 6e-4 0.0452 
Tempe WFon aw, see San 
77 413 595.8 33.1 0.569 235.9 251 


II]. MULTISTAGE FLASH VAPORIZATION 


Our main objectives were to recover as much as liquid 
condensate as possible while keeping its RVP at 10 psia for 
the summer case. However, heat duties of steam heaters and 
shell and tube heat exchangers (so as their surface areas) 
should be minimized to reduce the operating and capital costs 
to the least possible amount. The liquid stream leaving the 
second separator (L3 stream in Fig. 1) is recycled to preheat 
the liquid stream leaving the first one (L1 stream in Fig. 1) 
aiming at lowering the heat duty of the following steam 
heater (E-102 in Fig. 1). A pressure drop of 5 psi was 
designated to both shell and tube sides. The heat duty of 
E-102 is adjusted aiming at bringing the RVP of the final 
liquid product to 10 psia. The pressure of the first stage is set 
to the flow line pressure i.e. 139 psia (in Khangiran case). 
The storage tank pressure is 13.92 psia. The rule of thumb 
proposed by Campbell (2004) suggests that an equal pressure 
ratio between the stages would be considered. Here, it gives a 
pressure of 44 psia for the second separator (V-102 in Fig. 1). 
However, the trend of the total heat duty (the sum of heat 
duties of E-100, E-101, and E-102 heat exchangers) and the 
recovered stabilized condensate shown in Figs. 2 and 3 


' Gas To Condensate ratio, STD m’.h" gas flow per STD m’.h" liquid 
flow 
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suggest that a lower pressure may be more economic. Hence, 
the second separator pressure was assumed 25 psia. In 
practice, a higher pressure should be selected to compensate 
for the pressure losses in the transport pipeline and to allow 
recovering the flashed gas in the second separator and 
sending it to the fuel gas system on demand. 


Stab. yay 
Cond, a 
1 Stab. Cond. 2 


_ — 


>< 


VLV-101 


— 
P To Storage 
E-100 Tank 


L3 
Recycled 
to E-101 
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Fig. 1. Multistage flash vaporization scheme with the recycle of hot 
stabilized condensate. 
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Fig. 2. Total heat duty versus second separator pressure in Khangiran case 
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Fig. 3. Stabilized condensate flow rate versus second separator pressure 
in Khangiran case 


The RVP and H,S concentration of stabilized condensate, 
the stabilized condensate to feed ratio, and the total heat duty 
of heat exchangers are listed in table III. The dimensions of 
separators could be considered as 4.5 ft in diameter and 15.75 
ft in length using horizontal vessels with length to diameter 
ratio of 3 and retention time of 12 minutes. 

In NGL3100 case, however, use of two flash separators 
doesn’t yield the required RVP of the condensate and cannot 
be considered as an option. We added one separator and 
investigated the flash vaporization scheme. Despite the high 
heat duty applied, the stabilized condensate recovery was 
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very low (3.5%) and the H»S content was much higher than 
that of fractionation. We can conclude that flash vaporization 
is not practical in NGL3100 case. 


TABLE III: STABILIZED CONDENSATE SPECIFICATIONS IN KHANGIRAN AND 


NGL3100 CASES 
Khangiran RVP (psia) HS (ppm) CFR’ (%) ~~ THD? (btu.h) 
Two Flash Tanks 10 4467 87 5.730 x10° 
Fractionation 10 6635 10 8.484 x10° 
NGL3100 
Three Flash Tanks 10 206 3.5 5.438x10° 
Fractionation 10 <l 30 5.639x10° 


IV. FRACTIONATION 


A distillation column with a reboiler and without any 
condenser (non-refluxed column) was used (Figs. 4 and 5). 
The pressure of top tray and reboiler are primarily considered 
as 130 and 135 psia in Khangiran case and 140 and 145 psia 
in NGL3100 case respectively. Boil-up ratio was adjusted 
aiming at having the bottom product RVP at 10 psia. The feed 
is preheated to an extent that the preheater, reboiler, and 
cooler (E-101 in Fig. 4) heat duties are lowered as much as 
possible in Khangiran case. It is notable that the stabilized 
condensate flow rate and the total heat duty variations have 
the same trend indicating that there’s a trade-off between the 
heating cost and the stabilized condensate profit. 

The criterion of preheating extent is different in NGL3100 
case. It is determined in such a way that the temperature of 
cooled stabilized condensate (DL2 in Fig. 5) reaches to 95 F 
(which is almost the stock tank temperature). Since the feed 
pressure of NGL3100 is high (413.4 psia) it passes through a 
Joule-Thomson valve (VLV-101 in Fig. 5) and reaches to 
130°F and 159.6 psia. 
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STABCOL 
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Fig. 4. Fractionation (Distillation) scheme with the recycle of hot 
stabilized condensate (bottom product) in Khangiran case 


For a definite separation, the reflux and boil-up ratios 
(consequently the total heat duty of condenser and reboiler 
and their respective surface area) decrease as the number of 
stages increases. However, increase in the number of stages 
results in higher capital costs due to higher heights of column 
and trays price. Therefore, there’s a trade-off between the 
number of stages and the total heat duty. On the other hand, 
the flow rate of stabilized condensate rises with respect to the 
number of stages. The changes in reboiler duty and flow rate 


? Condensate To Feed Ratio (Ibmol cond./Ibmol feed) 
> Total Heat Duty 
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of stabilized condensate in NGL3100 case are depicted in 
figures 6 and 7. It can be seen that the addition of stages 
further than 26 has small effect on the reboiler duty and the 
condensate flow rate. Hence, this number may be primarily 
considered as the optimized one. If the overall tray efficiency 
is considered as 0.85, the actual number of stages will be 31 
which equals to 30 trays. The valve tray type with one pass, 
flooding fractional approach of 0.8, and tray spacing of 1.5 ft 
were chosen. The calculated values for diameter and column 
pressure drop were 3.07 ft and 1.57 psi respectively. A 
similar procedure in Khangiran case indicates 10 theoretical 
and 13 actual stages. The column diameter and pressure drop 
were calculated as 2.17 ft and 0.76 psi respectively (at 1 ft 
tray spacing). Considering the diameter of column (2.17 ft), 
choosing a packed column would be better. If we use | in. (25 
mm) pall ring packing with approximate HETP* of 1.48 ft 
(0.45 m), the column diameter and pressure are calculated as 
2.02 ft and 0.04 psi respectively. The properties of product 
and the total heat duty of reboiler and heat exchangers are 
listed in table ITI. 
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Fig. 5. Fractionation (Distillation) scheme with the recycle of hot 
stabilized condensate (bottom product) in NGL3100 case 
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Fig. 6. Reboiler heat duty versus number of stages of distillation column 
in NGL3100 case 
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g. 7. Stabilized condensate flow rate versus number of stages of 
distillation column in NGL3100 case 
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More rigorous designs should take economic evaluations 
into account. The capital cost of column is related to the 
vessel, trays or packings, internals, and the reboiler (a kettle 
type shell and tube heat exchanger). This cost may be 
annualized via dividing by payback period in years. The 
operating costs account for the energy required for generation 
of high pressure steam passing through the reboiler tubes to 
partially vaporize the bottom liquid. Total annual cost (TAC) 
may be then calculated from (1). The number of trays which 
yields the minimum TAC may be considered as the optimized 
number of trays. 


TAC = capital cost 


, + energy cost (1) 
payback period 


The f.0.b. purchase cost of the distillation column can be 
calculated by (2) to (11) [5]. 


C,=C,+C, (2) 
60a ee On en (3) 


Cy =exp(7.0374+0.18255[In(W)]+0.02297[ In(w)? ]) (4) 


W =n(D, +t, )(L+0.8D,)t,.p (5) 

Ce = 285.1(D) (Dy (6) 
Cea IN gcP Ba Pace (7) 

we ——. Ne<20 (8) 

Ht N< 20 (9) 

Fy =1.401+0.0724 D, (316S.S.) (10) 
Cyr = 369 exp(0.1739 D,) (11) 


where, Cp is the distillation column cost ($), Cs is the shell 
cost ($), Cr is the tray cost ($), Fu is the material factor which 
is 2.1 for 316 stainless steel (S.S.), Cy is the cost of empty 
vessel including nozzles, manholes, and supports ($), Cpr is 
the platform and ladder cost ($),W is the weight of shell and 
two heads (lb), D; is the inside diameter (in), t, is the shell 
thickness (in), L is the column height (in), p is the 
construction material density (Ib.in*), Ny is the number of 
trays, Fr; accounts for the tray type which is 1.18 for valve 
trays, Fry is the correction factor for the material of 
construction, and Cpr is the base cost of trays. 

The f.o.b. purchase cost of the reboiler which is a kettle 
type shell and tube heat exchanger is calculated by (12) to (15) 


[S]. 


C, = Fy.Fy-F,.C, (12) 


P py 
F,=0.8510+0.1292 2}, 0.0198 2 | (13) 
600 600 


b 

A 

F, =a+| — 
ee (=) 


C= exp(I 1.967 —0.8709[In(A)] +0.09005[In(4)] (15) 


(14) 


where, Fp is the pressure factor, P is the shell side pressure 
(psig), Fy is the material factor, a and b are 2.70 and 0.07 for 
S.S./S.S. material, A is the surface area (ft’), and Cx is the 
base cost ($). Energy cost may be considered as 6 $.GJ". 
The aforementioned equipment cost formulas are indexed to 
the year 2000 and should be updated to the year 2011. It can 
be done by multiplying the base cost by the cost index ratio as 
stated in (16). We used Chemical Engineering (CE) index 
which is 394 and 585.7 for the years 2000 and 2011 
respectively. 


(16) 


I 
Current cost = base cost [i= 
2000 


The diagram of TAC with respect to the number of trays is 
shown in Fig.8. It implies that the number of trays at which 
the minimum TAC occurs is 9 (or 10 stages). 
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Fig. 8. Total annual cost versus number of trays of distillation column in 
NGL3100 case 


Since the stabilized condensate flow rate increases with 
respect to the increase in number of stages, it is better to 
choose the optimized number of trays by maximizing the 
annual profit stated in (17). Annual sale is calculated by 
multiplication of gas condensate price (in $.bbI') by gas 
condensate flow rate (in bbl.y”'). The stabilized condensate 
price was taken as 100 $.bbI’'. For simplicity, we calculated 
the annual profit obtained by the addition of each tray (profit 
per tray) instead, by assuming the number of two trays as a 
base. Since the duty and area of the preheater and the 
associated equipments do not change significantly with 
respect to the addition of column trays, the aforementioned 
trick would be rational. Fig. 9 shows the annual profit 
obtained by the addition of each tray. It’s obvious that profit 
eared by using more than 9 up to about 30 trays is 
considerable. Hence, the number of 30 trays (or 31 stages) 
which was suggested earlier sounds reasonable. 


Annual profit=Annual stabilized condensate sale-TAC (17) 
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Fig. 9. Annual profit per tray versus number of trays of distillation 
column in NGL3100 case 


V. PROCESS SELECTION 


We saw that the multistage flash vaporization is a practical 
method to stabilize the condensate in Khangiran case, as it is 
already utilized in the aforementioned gas processing plant. 
However, in NGL3100 case, it doesn’t sound promising 
considering the low condensate recovery and high total heat 
duty. Capital costs of the process equipments aside, 
multistage flash vaporization is more suitable than the 
fractionation in terms of total heat duty and the flow rate of 
stabilized condensate (or condensate to feed ratio) in 
Khangiran case. However, in NGL3100 case, it is clear from 
table III that fractionation is far better than flashing since it 
provides more condensate recovery, less heat duty, and H2S 
content within the standard range. It seems that RVP of the 
feed has a great effect on the performance of stabilization unit 
which, in these cases, were 155.7 psia and 235.9 psia. We 
should take this fact into account that distillation column is 
further capable of handling variations in the feed composition 
(which affects its RVP). H2S content is high in both methods 
in Khangiran case and must be reduced to standard values 
(e.g. less than 50 ppm). 


VI. CONCLUSION 


In the present study, two methods of multistage flash 
vaporization and  non-refluxed fractionation were 
investigated to stabilize gas condensate by reducing its RVP 
to 10 psia in summer. Gas condensates of two gas processing 
plants were studied. Results show that depending on the feed 
properties (especially its RVP) one or both of the 


aforementioned techniques is/are applicable and/or economic. 


In Khangiran case (RVP and GCR of the feed were 155.7 psia 
and 179 respectively), multistage flashing provides less total 
heat duty, more condensate recovery and less HS content 
which makes it more interesting. In NGL3100 case (RVP and 
GCR of the feed were 235.9 psia and 251 respectively), 
fractionation yields less total heat duty and much higher 
condensate recovery. It is also capable of lowering H2S 
concentration to less than 1 ppm. Therefore, fractionation is 
highly preferred. Calculations made to determine the 
specifications of distillation columns and flash separators as 
follows: a non-refluxed trayed distillation column with 31 
stages (30 trays), 3.07 ft in diameter, and 54 ft in height in 
NGL3100 case; and a similar type of column with 13 stages 
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(12 trays), 2.17 ft in diameter, and 14.4 ft in height or a 
packed column with the same number of stages but 2.02 ft in 
diameter and 21.3 ft in height in Khangiran case; and two 
cylindrical horizontal separators with dimensions of 4.5 ft in 
diameter and 15.75 ft in length in Khangiran case. Finally, 
detailed cost benefit analysis together with safety 
considerations and environmental concerns would determine 
the final selection of the right technique. 
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